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Wc introduce the concept of back-reaction in relativistic cosmological modeling. Roughly 
speaking, this can be thought of as the difference between the large-scale behaviour of 
an inhomogeneous cosmological solution of Einstein's equations, and a homogeneous and 
isotropic solution that is a best-fit to either the average of obscrvablcs or dynamics in 
the inhomogeneous solution. This is sometimes paraphrased as 'the effect that structure 
has of the large-scale evolution of the universe'. Various different approaches have been 
taken in the literature in order to try and understand back-reaction in cosmology. We 
provide a brief and critical summary of some of them, highlighting recent progress that 
has been made in each case. 
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1. What is Back- Reaction? 

The term 'back-reaction' is often used in cosmology to mean 'the effect that struc- 
ture has on the large-scale evolution of the universe, and observations made within 
it'. Implicit within this statement are a number of fundamental problems that have 
yet to be fully understood. These include: 

1. What is meant by the large-scale expansion of space in an inhomogeneous 
universe, and how should it be calculated? 

2. How should we link the large-scale expansion of an inhomogeneous space- 
time with the observations made within it? 

3. How can we create relativistic cosmological models sophisticated enough to 
investigate these problems? 

Let us now briefly consider each of these points, before moving on to discuss recent 
attempts to understand them. 

Point 1 above alludes to the fact that in relativistic theories what we mean by 
the spatial separation of any two astrophysical objects depends on how we choose 
to foliate the universe with hyper-surfaces of constant time. In a spatially homoge- 
neous universe, or a universe with an irrotational matter content, natural-looking 
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choices might present themselves. In general, however, we should be free to make 
any number of choices. This then presents a problem: If the distance between any 
two astrophysical objects is in general foliation dependent, and we have no preferred 
foliation, then how should we go about defining the rate of change of distance be- 
tween objects, and hence the expansion of the universe? In the end, the answer to 
this question will depend on exactly what one is trying to achieve, and is compli- 
cated considerably by the fact that in cosmology one is often interested in non-local 
averages (a notoriously difficult concept to define in general relativity). Below we 
will consider several different cases of interest. 

Point 2 is a subsequent problem that needs to be addressed, once a concept of 
'large-scale expansion' exists that one is prepared to consider. It is not in general 
the case that observations made in an inhomogeneous geometry will have a straight- 
forward correspondence with the observables that one would measure in a spatially 
homogeneous and isotropic universe with the same rate of expansion on large scales. 
That is, even if one succeeds in finding a good description for the large-scale expan- 
sion of the universe, then one still needs to do further work in order to relate this to 
observations made in the underlying inhomogeneous space-time. Once again, this 
is complicated considerably by the fact that we are often interested in the average 
of observables. This is in general a highly non-trivial problem, and below we will 
review some recent progress towards understanding it. 

Finally, point 3 is related to the fact that in order to test proposed solutions to 
the problems posed in points 1 and 2 it is of considerable interest to have cosmolog- 
ical models that are sophisticated enough to allow at least some of the interesting 
behavior that we expect in general. This is an extremely difficult problem. Although 
many inhomogeneous cosmological solutions to Einstein's equations are known,E 
most of these solutions are restricted either because they are required to exhibit a 
high degree of symmetry or because they are algebraically special. Constructions 
such as the "Swiss cheese" models allow some potential progress to be made, but 
are themselves severely restricted by the boundary conditions at the edge of each 
"hole" . New approaches are required to make further progress in this area, and, 
once again, we will discuss some recent progress below. 

In Section [2] wc consider approaches based on averaging over a set of prescribed 
spatial hyper-surfaces. In Section [3] we consider approaches based on averaging in 
four dimensions. Section [4] contains a discussion of some models that may be of use 
for studying averaging, and in Section [5] we provide a few closing comments. 

2. Spatial Averaging Approach 

One way to proceed with the study of back-reaction is to consider the expansion 
of regions of space in a given foliation. The equations that govern this expansion 
can then be found, and compared to the Fricdmann equations. This often leads one 
to consider the volume-weighted average of quantities such as energy density and 
pressure. The equations that result are therefore often referred to as the 'averaged 
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field equations'. 

While simple, this approach has a number of obvious drawbacks. Firstly, it is 
manifestly not foliation invariant. Secondly, there is a freedom in how one chooses 
to specify that two spatial volumes at different times are the same region. And 
thirdly, the averaging of quantities over the spatial volume being considered is often 
only well defined for scalars. One can specify choices for the first and second of 
these that may initially appear natural, but that could in the end lead one to 
consider hyper-surfaces in the inhomogeneous space-time that become arbitrarily, 
and increasingly, distorted. The third of these problems is of more fundamental 
difficulty, as tensors cannot in general be compared at different points. Nevertheless, 
this approach provides a useful framework to investigate, and can be shown to give 
a straightforward correspondance to the average of observables in some situations. 



2.1. Buchert's Equations 

The most well studied set of averaged equations that result from this approach 
are those found by Buchert after averaging the Hamiltonian, Raychaudhuri and 
conservation equations. 2 

^ = ^ n{p ) v -1(^r) v -Iq v (1) 

3^ = -4<kG n ( P )v + Qv (2) 
a-D 

d t (p)v + S—(p)v = 0, (3) 
an 

where a-D and (^R) are the "scale factor" and average Ricci curvature of the region 
of space T> being considered, angular brackets denote a volume average throughout 
that region, and Qd is the back-reaction term that quantifies differences from the 
Friedmann equations that one might otherwise construct from these quantities. 
These are defined as 

\S„<PXy/Wg(to,X*)) 



J^X^^W^) 



where and a are the expansion and volume-preserving shear of the set of curves 
orthogonal to the hyper-surfaces containing T>, and t and X 1 are the proper time 
measured along this set of curves and the spatial coordinates in the hyper-surfaces 
of constant t, respectively. The quantity to is the value of t on some reference hyper- 
surface (usually taken to be the one that contains us at present). 
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One may note that equations do not form a closed set. Extra information 

is therefore required, which can be given by specifying Qx> = Qv{t)- Presumably 
this requires either extra equations, or some knowledge of the inhomogeneous space- 
time being averaged. As previously stated, one may also note that the averaging 
procedure given here by the angular brackets is foliation dependent and only appli- 
cable to scalars (this is particularly problematic for the term (<r 2 )x> in equation M, 
as the evolution equation for a 2 will contain tensors). Finally, while the expansion 
of the spatial domain T> may not itself be directly observable, we will explain below 
that in some cases it can be linked to observables. 



2.2. Links to Observables 

The term 'observables' can cover a wide array of different possibilities in cosmology. 
Here we will mainly be concerned with the luminosity distance-redshift relation. 
This is itself a direct observable of considerable interest for the interpretation of, 
for example, supernova observations. Beyond this, it is also often required in the 
interpretation of other observables as it is very often the case that one needs to 
transform from "redshift space" to some concept of position space (i.e. the position 
of astrophysical objects on some spatial hyper-surface). 

The usual method for calculating luminosity distances in an inhomogeneous 
space-time is to first find the angular diameter distance to the emitting object as 
a function of some affine parameter, measuring distance along past-directed null 
geodesies. This can be achieved by integrating the Sachs optical equations™ In 
these equations the Ricci curvature of the space-time sources the evolution of the 
expansion of the past-directed null geodesies, and the Weyl curvature sources the 
evolution of their volume-preserving shear (which itself acts as a source for their 
expansion). The angular-diameter distance can then be straightforwardly related 
to the luminosity distance,^ and the redshift can be calculated as a function of the 
affine distance (once the world-lines of the objects emitting the radiation have been 
specified). This then provides the luminosity distance as a function of redshift at 
all points on an observer's past-light cone, provided that geometric optics remains 
a good approximation, and that the light emitted from the distant object is not 
obscured by some intermediate matter before it reaches the observer. 

Although the method outlined above is, in general, a complicated problem in- 
volving a number of subtleties, it was recently shown by RasanerPthat progress can 
be made in space-times that display statistical homogeneity and isotropy on large 
scales. In this case one can estimate the average luminosity distance as a function 
of the average redshift that an observer in such a space-time may expect to recon- 
struct from observations made over cosmologically interesting distances. Assuming 
that the matter content is irrotational, that the shear in the null trajectories can 
safely be assumed to be small, that structures evolve slowly, and that hyper-surfaces 
of constant proper time can also be taken to be the same hyper-surfaces that display 
statistical homogeneity and isotropy, Rasanen made a convincing case that the av- 
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erage luminosity distance-redshift relation in the inhomogeneous space-time should 
be well approximated by observables calculated in a homogeneous and isotropic 
model with a scale factor that evolves according to equations ([TJ)- 

An alternative approach to this problem was taken by Clarkson and UmehP 
These authors considered expressions for measures of distance expanded as a power 
series in redshift, as derived for general space-times by Kristian and Sachs.^ They 
then performed a decomposition into spherical harmonics, and constructed the fol- 
lowing deceleration parameter, based on an analogy between the monopole of this 
expansion and the corresponding relations in a Friedmann universe: 



(7) 

o 

where H = 0/3 is the isotropic part of the Hubble rate, and subscript '0' denotes a 
quantity evaluated at z = 0. Using this expression they could consider the average 
deceleration within cither a region of space, or a region of space-time. However, for 
matter obeying the strong-energy condition it can be seen from equation Q that 
the average of qo will always be non-negative, and so the space-time (according to 
this measure) will always be inferred to be decelerating (in the absence of A). This 
is in contrast to the averaged evolutions possible from equations Q-Q, and at first 
glance would appear to contradict the results of Rasanen described above. 

In fact, there is no contradiction between these two sets of results P That is, the 
observable calculated by Clarkson and Umeh should be expected to be a good ap- 
proximation to the deceleration that one would infer from observations made within 
a small region around an observer. This measure is closely related to the acceler- 
ation of space within that region, as specified by Einstein's equations (as long as 
shear is small), and not by equations ([l])-([3]). The observational measures consid- 
ered by Rasanen, however, are only expected to approach the evolution described 
by equations ([l])-([3]) when the distances over which observations are made are much 
larger than the homogeneity scale of the space-time under consideration. This is, 
of course, the regime in which cosmological observations are usually made. Using 
example space-times it has been explicitly demonstrated that it is entirely possible 
for a set of observers in a given region of the universe to infer deceleration from 
Clarkson and Umeh's measure, while inferring acceleration from Buchert's mea- 
sureP This clearly demonstrates that the acceleration inferred from cosmological 
observations does not have to be closely related to the local acceleration of space 
itself. It also demonstrates that quantities that are uniquely defined in an exactly 
homogeneous and isotropic universe (such as go) can bifurcate into multiple differ- 
ent quantities in space-times that are only statistically homogeneous and isotropic, 
and that in general these new quantities can take very different values from each 
other. One must therefore proceed with care. 



5o = 
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3. Space-Time Averaging Approach 

An alternative approach to considering the volume weighted average of quanti- 
ties within 3-dimensional spatial regions is to instead consider averaging geometric 
quantities within 4-dimensional regions of space-time. Such a process is in general 
difficult to define in a covariant way, and so far has required the application of 
bi-local operators. These allow tensors to be compared at different points by trans- 
porting them along prescribed sets of curves. This then leads to the problems of how 
the curves in question should be prescribed, and exactly which transport method 
should be used. Various proposals exist as to the best way to address these issuesP 
While complicated, the idea of averaging quantities in 4-dimensional regions of 
space-time inherently avoids any foliation dependence. These approaches are also 
often aimed at averaging tensors directly, rather than just scalars. This has obvi- 
ous advantages for gravitational theories constructed from tensors, such as general 
relativity. 



3.1. Zalaletdinov's Equations 

Probably the most well known attempt at averaging in space-time, and constructing 
a set of effective field equations that the averages should obey, is that of Zalalet- 
dinov!^ The first step in this approach is to construct the following average for a 
tensor Pg'y: 

(pf.rix)) = JL jf y^V^; : >vV(^V'>,*') (8) 

where primed coordinates are those used in the 4-dimensional region S, which is 
the averaging domain associated with the point x. The quantities A a ^, are the 
bi-local operators, which are functions of both x and x' , and the quantity Ve = 
Js ^/—g'd 4 x' is the volume of S. Each point, x 7 is expected to have associated with 
it its own averaging domain, S, which is related to other averaging domains by 
being transported around the manifold. 

By applying this averaging technique to the connection, and by using some 
"splitting rules" , Zalaletdinov is able to use Einstein's equations to derive a set of 
field equations that the averaged connection must obey. The are called the Macro- 
scopic Field Equations, and are writterP^ 

gfi'M^p - \s%g^M^ = 8nGT% - (z*^ - \s%Q^ gT , (9) 

where bars denote averaged quantities, and M 1 p = M a ^ a/3 and Q^ v = Z a ia/a and 

^ = 2Z V^i' where 

M" vafl = d a (v\ p ) - d p (T^ a ) + (r^ a )(v% p ) (r%)(r%j (10) 
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and where underlined indices are not included in symmetrization operations. The 
tensor Z a ^ ua is known as the 2-point correlation tensor, and obeys its own alge- 
braic and differential constraints.^ 

The Macroscopic Field Equations (|9| can be used to describe the behavior of a 
particular inhomogeneous space-time after averaging has been performed, but they 
can also be used as a set of field equations to which one can look for solutions 
directly. This latter approach has so far been taken in the cases of macroscopic 
geometries, g M „, that are spatially homogeneous and isotropic}iiJli2l and geometries 
that are spherically symmetric and static.^ This work has allowed some possible 
behaviors of averaged space-times to be found without specifying the underlying 
microscopic geometry. However, it has also so far required a number of assumptions 
to be made about the correlations that are present. These include the vanishing of 
the three-point and four-point correlation tensors, and the vanishing of the 'electric' 
part of the 2-point correlation tensor.^ The particular situations in which these as- 
sumptions are valid remains to be determined, as is also the case for the assumptions 
that go into the derivation of the Macroscopic Field Equations ([9]). Nevertheless, 
this is an interesting approach that deserves further study. 



3.2. Links to Observables 

Under the assumption that the macroscopic geometry is spatially homogeneous and 
isotropic (that is, after the averaging procedure has been applied, and the "aver- 
aged" geometry displays these symmetries), then Coley, Pelavas and Zalaletdinov 
find the following to be a solution of the Macroscopic Field Equations Q: 11 

r 2 (d0 2 + sin 2 Odcf) 2 )] , 



g^dx^dx" = -dt 2 + a 2 (t) 



1 — k g r 2 



(12) 



where k g is a constant, and where a(t) and p obeys the Fricdmann-like equations 

a 2 8nG k d 

^ = ^r p -^ (13) 

p + '3-(p + p)=0, (14) 
a 

where kd is a constant (not necessarily equal to kg), and where p and p are the 
macroscopic energy density and pressure (obtained after averaging the right-hand 
side of Einstein's equations). 



Superficially, the geometry given in equations ( 12 1-( 14 ) looks a lot like the spa- 
tially homogeneous and isotropic solutions of Einstein's equations in the presence 
of a perfect fluid. There is, however, a very significant difference: The spatial cur- 
vature constant that appears in the macroscopic geometry, k g , is not in general the 
same as the term that looks like spatial curvature in the Friedmann-likc equation 



(13) (i.e. the one that contains fc^). That is, spatially curvature can take different 
values depending on the situation being considered. If one measured the angles at 
the corners of triangle, and determined the curvature of space in this way, then 
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this would give a different result to that which would be obtained by measuring 
the recessional velocities of astrophysical objects and inferring the spatial curvature 



through the dynamical (Friedmann-like) equation (13 1. This behavior is impossible 
within the spatially homogeneous and isotropic solutions of Einstein's equations, 
and so could provide some potentially observable phenomena that could be used 
to test this approach. The difference between k g and kd is determined by terms 
that appear in the correlation tensor, Z a ^ l V(y , and so by attempting to determine 
the difference between k g and kd observationally we could attempt to constrain 
Z a gJ l ua , and hence some of the possible effects of averaging. 

A first step towards investigating this possibility has recently been takenP^l 
The authors of this work assume that average observables are determined by null 



trajectories in the average geometry, as specified in equation ( 12 ), and that redshifts 
are represented by the average scale factor, a(t). They then find that luminosity 
distances are given by the following equation: 

■ n= (1 + *) . ( f 1 VW^\da \ 

L[Z> #oVTO \Jjh VV kd a 2 + Q A a 4 + n m a) ' 1 ' 

where the matter content of the macroscopic space-time has been assumed to be 
well approximated by non-interacting dust and A, where Hq — a/a\ z =0: an d where 
the fL are defined as 



o _ k a o - kd o _ 8nGp„ h0 _ 87rGp A 

il k g — 72772' "Aid— 72772' "m — ^772 , "A - „ „ 2 ' \ W ) 

a n "o-^o «" 3 a 



where p m ,o is the present energy density in dust, and 87rGpA = A is the effective 



energy density in A. The expression for luminosity distance given in equation ( 15 ) 
can now be used to interpret cosmological observations, and to obtain constraints 
on the fli. 

Using data from the Hubble Space Telescope (HST)f^the Wilkinson Microwave 
Anisotropy Probe (WMAP)J^ observations of the Baryon Acoustic Oscillations 
(BAOs)p^and the Union2^and SDSS^I supernova data sets, the parameters flk d , 
flk and r?A were constrained to take the values given in Table [l] belowP The ad- 
ditional freedom of allowing Q,k d ^ f2& in this analysis means that the CMB+i?o 
is now no longer sufficient to constrain the spatial curvature of the universe signifi- 
cantly. Observations of the CMB+Hq alone are also no longer sufficient to require 
A ^ 0. This simple extra degree of freedom therefore undermines two of the most 
important results of modern observational cosmology. By adding further data sets 
the constraints on £lk d and flk g are improved, but still remain much weaker than in 
the standard Friedmann models that satisfy Einstein's equations. Even so, however, 
it was still found that the results of using all available observables were sufficient to 
require f^A 7^ to high confidence, and that a spatially flat universe was consistent 
with observations. Finally, although the combination of some data sets excluded 
the possibility f2/c d = £lk g at the 95% confidence level, it was found that the special 
case Qk d = fl^ was compatible with most combinations of these data sets. 
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Table 1. Constraints on f^fe^, &k g snd from data sets outlined in the text. 



Data Sets 


«*«, 






CMB 
CMB+HST 
SNIa (Union2) 
SNIa (SDSS) 
CMB+HST+SNIa(Union2) 
CMB+HST+SNIa(SDSS) 
CMB+HST+SNIa(Union2)+BAO 
CMB+HST+SNIa(SDSS)+BAO 


-0.053t£i£ 
+0.036ig;gg 
+0.0121°;™ 
+0.233irj? 

+o.oi4±8;gi? 

+0.054ir 2 2 ° 
-0.004±g;gii 
+0.026lHiI 


-0.036+^^ 
+0.185l°if 5 
-0.369™ 

-o.i73i°;f 2 7 

+0.055tr f 2 

in •311+0-100 
-|-u.ciii_ 101 

-o n33+ - 070 

u.uoo_ 069 

4-n loq+0.072 


+0 ■525+ - 417 
+0.564+.°;^ 

+0.902™ 
+0.641+ i° 5 

in fiq^+0.080 
-r-u.oao_ 0g2 

_i_n 4 q fi +0.087 

+o.75B±8;g?g 

+0.522ir 7°3 



In general one might also consider the possibility of not just allowing ttk d and fit 
to be different, but also allowing them to functions of scale. Such a result might arise, 
for example, from performing averaging over domains of different sizes, a process 
which is implicitly carried out when consider different cosmological observables. 
Such a possibility allows for considerable extra freedomFS 

4. Constructing Inhomogeneous Models 

We have so far considered attempts to describe the large-scale behavior of the uni- 
verse by averaging over regions of space or space-time. In the end, the particular 
approach that one should use when performing this type of operation should prob- 
ably be guided by the phenomena that one is trying to create a model to interpret. 
Different observable phenomena may require different approaches, and so one needs 
to know the limits of any particular approach, as well as the situations in which it 
reliably reproduces the required results. For this it is useful to have inhomogeneous 
cosmological models that are of sufficient generality to allow some of the interesting 
behavior that is expected in general. Such models can then be used to test ideas 
about averaging, back-reaction, and the large-scale evolution of space. 

Unfortunately it is extremely difficult to construct such models. This does not 
mean that there is an absence of any interesting behavior to study, only that we need 
to become more sophisticated in our model building to quantify and constrain the 
different possibilities in a reliable way. Some of the principal difficulties involved 
with this are how to model over-dense regions of the universe without having to 
deal with the rapid formation of singularities, how to introduce structure into the 
universe without assuming a Friedmann background or matching onto a Friedmann 
model at a boundary, and how to allow structure to form on different scales without 
assuming linearity in the field equations. For further discussion of inhomogeneous 
cosmological solutions the reader is referred to the contribution to these proceedings 
by KrasihskipSI and to the comprehensive textsEE] 

It is currently almost beyond hope to construct a model that allows for all of the 
possibilities discussed above, while simultaneously maintaining sufficiently general- 
ity to model realistic distributions of matter. We are therefore forced to investigate 
toy models that we hope may reflect some of the features of the real universe, even 
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if they are not realistic in every way. Once toy models have been constructed we 
can then consider the averaging problem by applying some of the methods discussed 
above to them, or by fitting or comparing them to Friedmann models directly. Their 
existence also makes more advanced models a more realistic proposition. It is for 
these reasons that it is of interest to consider simple n-body solutions of Einstein's 
equations. Such solutions, if they can be found, will allow over-dense regions to be 
studied without rapid collapse occurring, and without recourse to the assumption 
of a Friedmann background or linearity in the gravitational field equations. This 
will be the subject of Section |4~T) 

4.1. A Lattice of Black Holes 

The simplest configuration of n bodies that one can imagine is a regularly arranged 
set of points. Although such a configuration limits the behaviors that are possible, 
it does allow for the most straightforward possible comparison to smoothed-out 
Fricdmann-like universes. That is, by 'zooming out' in order to consider large num- 
bers of points, and by performing some kind of coarse graining or smoothing, one 
could easily imagine such a situation looking more and more like a spatially homo- 
geneous and isotropic universe, which could then be compared to the Friedmann 
solutions of Einstein's equations. Regularity of the distribution also provides a lim- 
ited number of preferred spatial planes and curves that can have their area and 
length compared to those of the Friedmann solutions. 

Here we will consider spatially closed universes. These are known to admit hyper- 
surfaces of time symmetry at the maximum of expansion of the space-time that allow 
the constraint equations to be solved in a particularly simple wayP^l The method 
that we will deploy to ensure that our massive bodies are regularly arranged is to 
tile the hyper-surface of maximum expansion with a number of regular polyhedra. 
A mass is then placed at the center of each polyhedron, which by symmetry will 
be an equal distance from each of its nearest neighbors. These polyhedra will be 
referred to in what follows as 'cells'. There are seven such tilings that are possible 
in three spatial dimensions, as listed in Table [2j Also displayed in this table are the 
Schlafli symbols of the polychora that these tilings constituted 



Table 2. Tilings of the 3-space of maximum expansion, and their scale in 
comparison to the homogeneous and isotropic Friedmann solutions. 



Lattice 


Cell 


Number of 


Ratio of scales in discrete 


Structure 


Shape 


Cells 


and Friedmann solutions^! 




Ball 


2 




{333} 


Tetrahedron 


5 


1.360 


{433} 


Cube 


8 


1.291 


{334} 


Tetrahedron 


16 


1.097 


{343} 


Octahedron 


24 


1.099 


{533} 


Dodecahedron 


120 


1.034 


{335} 


Tetrahedron 


600 


1.002 
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Once the arrangement of masses has been chosen, the geometry of the hyper- 
surface of maximum-of-expansion can be found. With the exception of the 2-cell, 
this has been done for each of the structures described above The 2-cell is special 
in that the time-symmetric geometry at the maximum of expansion of this struc- 
ture is simply a slice through the global Schwarzschild solution. The geometry of 
the full space-time is therefore already known exactly in this case, and is not of 
cosmological interest here (by including a non-zero A, however, other structures are 
also possible^). An illustration of the geometry at the maximum of expansion in 
the case of the 8-cell and the 120-cell is given in Figure [l] below. Each of the illus- 
trations here corresponds to a single 2-dimensional slice through the 3-dimensional 
geometry. In the case of the 8-cell this slice contains 6 masses, while in the 120-cell 
it contains many more (although not all 120). The geometry of the 3-space of maxi- 
mum of expansion in each case is conformally related to the geometry of a 3-sphere, 
with a scale factor that is a function of position. The distance from the origin in 
the illustrations in Figure [T] is proportional to this scale factor, and it can be seen 
that as the number of the masses in the lattice is increased, the bulk of the space 
approaches homogeneity. It is only in the vicinity of the masses themselves that 
inhomogeneities exist (as depicted by the tube- like structures). 




(a) A slice through the 8-cell. (b) A slice through the 120-cell. 



Fig. 1. Graphic illustrations of the geometry of space at the maximum of expansion. 

Once we have the geometry of the hyper-surface of maximum of expansion, we 
can take a measure of the scale of the solution, and compare this to the scale of 
a spatially closed Friedmann universe that contains the same amount of "proper 
mass^'. The Friedmann solutions will, of course, have this mass evenly distributed 
throughout space, and so by comparing to the scale of the inhomogeneous geometry 
we can obtain a measure of back-reaction. For the choice of scale in the inhomoge- 
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neous space on could choose a number ol different measures. Here we consider the 
proper length of the edge of a cell. This corresponds to the scale of curvature for 
the sphere that appears to emerge when the number of masses becomes large (as 
can be seen from the illustration in Figure lb). The difference in scale in each case 
is given in the last column of Table [2j It can be seen that the broad trend is for the 
scale of the homogeneous and inhomogeneous space-times to approach each other 
as the number of cells becomes large. However, for only a small number of cells 
(~ 5 to 24) the difference in scale can be of order 10%. In any case, this method 
provides an exact quantification of back-reaction, and provides an arena for testing 
formalisms designed for more general configurations of energy and momentum. 

A numerical evolution of the 8-cell has now been performed,^ and other methods 
have also been used to address problem of understand the evolution of this type of 
structurePtH 



5. Discussion and Outlook 

Various approaches to back-reaction and averaging already exist in the literature, 
but much work remains to be done if we are to fully understand their observational 
consequences in the real universe. Motivation for taking these problems seriously 
comes from the apparent necessity of including dark energy when we interpret obser- 
vations within a linearly perturbed Friedmann model, as well as the requirement to 
understand all possible sources of error and uncertainty in precision cosmology. To 
fully address this problem it is likely that we will need to develop more sophisticated 
models of inhomogeneous space-times, as well as developing a more sophisticated 
understanding of averaging in general relativity. Research in this area should be 
considered exceptionally timely, with large amounts of resources currently being 
invested into observational probes designed to improve our understanding of dark 
energy, and the universe around us. 
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